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Abstract 

We exhibit a three parameter infinite family of quadratic recurrence relations 
inspired by the well known Somos sequences. For one infinite subfamily we prove 
that the recurrence generates an infinite sequence of integers by showing that the 
same sequence is generated by a linear recurrence (with suitable initial conditions). 
We also give conjectured relations among the three parameters so that the quadratic 
recurrences generate sequences of integers. 

1 Introduction 

It is always the case that quadratic recurrence relations generate sequences of rational 
numbers, given that the first finitely many terms are 1. That there exist quadratic recur- 
rences which generate only integers is very surprising. Michael Somos, through investi- 
gations of elliptic theta functions, discovered one such quadratic recurrence which began 
a new wave of research into this phenomenon [1]. After Somos' discovery [6J that the 
recurrence 

S{n) = (^S{n - l)S{n - 5) + S{n - 2)S{n - 4) + S{n - Sf^ /S{n - 6) 

generated integers when given initial conditions S{i) = 1 for < i < 5, a surge of 
research began on related sequences. Most importantly for the goals of our problem was 
the sequence discovered by Dana Scott defined by the recurrence 

D{n) = (^D{n - l)D{n - 3) + D{n - 2)^ /D{n - 4) (1) 

with initial conditions D{i) = 1 for < i < 3. This recurrence differs from Somos' original 
recurrence as it mixes quadratic and linear terms. However, it also defines a sequence 
consisting only of integers. 
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The above two sequences possess an even stronger property than integrahty. If the 
initial conditions are given as S{m) = Xm (resp. D{m) = Xm) for < m < 5 (resp. < 
m < 3), these recurrences produce Laurent polynomials. The set of Laurent polynomials 
are defined in [7] as 

M [t, t'l] = I ^ aif I i G Z, 7^ for finitely many i E z| 

We will define Laurent polynomials in more than one variable in the same manner as 
above. That is, the set of Laurent polynomials in the variables {xj}"^^ is defined to be 

K[{x.}r=i,{xri}r=i] = (2) 

{ ( n "^i^ ) I ^' £ Z, ttj 7^ for finitely many i G z| 

V i=i ^ 

This Laurent property is discussed in Fomin and Zelevinsky's article and has been 
used to extract combinatorial information from sequences. 

The work done by Dana Scott was important to our discovery since it gave us the idea 
to probe for sequences defined by a mixture of quadratic and linear terms, rather than 
just the Somos-like quadratic recurrences. Additionally, work done by Reid Barton [1] on 
a similar sequence resulted in a combinatorial interpretation, which led us to believe that 
a similar interpretation of our sequences exists. 



2 Statement of Problem 

While looking at the family of quadratic recurrences given by 

a{n)a{n - k) = a{n - \)a{n - A; + 1) + a(n - (A; - 1) /2) + ain - (A; + 1) /2) (3) 

with initial conditions aim) = 1 for m < k we noticed that they seemed to produce a 
sequence of integers for any odd value of k we chose. By substituting k = 2K + 1 we can 
transform the recurrence so it is valid for any positive integer value of K. This defines 
the following recurrence relation 

a{n)a{n - {2K + 1)) = a{n - l)a{n - 2K) + a{n - K) + a{n - K - 1) (4) 

now with initial conditions a(m) = 1 for m < 2K + 1. 
For instance, when K = 1 the recurrence becomes 

a{n)a{n — 3) = a{n — l)a{n — 2) + a{n — 1) + a{n — 2) 

generating the sequence {a{n)} = {1, 1, 1, 3, 7, 31, 85...}. 

We have proven the integrality assertion by showing that the sequence produced by the 
quadratic recurrence also satisfies a linear recurrence with integer coefficients and integer 
initial conditions. It has also been shown that when the initial conditions are given as 
formal variables the recurrences produce Laurent polynomials (see [5]). Our eventual 
goal is to find a combinatorial interpretation (e.g. counting perfect matchings of certain 
graphs) . 
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3 Proof of Integrality 



In order to prove that the quadratic recurrence produces integers we show that the se- 
quence also satisfies a hnear recurrence with integer coefficients. 

Theorem 1. The recurrence 

a{n)a{n - {2K + 1)) = a{n - l)a{n - 2K) + a{n - K) + a{n - K - 1) 
with initial conditions a{m) = 1 form < 2K + 1 generates an infinite sequence of integers. 

Lemma 1. The initial 6K + 1 terms of the sequence given by the quadratic recurrence 
relation Ol) satisfy the following conditions: 



a{2K + i 
alsK + i 
a{AK + i 
a{5K + i 



1 l<i<2K +1 

-l + 2i 2<i<K+l 

1 + 2K -2i + 2i'^ 2<i<K+l 

-3-8K + 2i- 2K^ + l2Ki + 2i'^ + AKH 2<i<K + 1 
3 + im - m + IQE"^ - IQKi + 8i^ 

+AK^ - AKH + 16Ki^ + AKH'^ 2<i<K + l 



Proof. Notice that a{i) = 1 for 1 < i < 2K + 1 by definition. Each of the other relations 
is proved independently by induction. All initial relations were originally conjectured and 
proved via a computer program written by Doron Zeilberger [8j. We will now show the 
proof for a{2K + i) = -1 + 2i with 2 < i < K + I. 

Base case i = 2: We need to verify that a{2K + 2) = 3. From the quadratic definition 
of the sequence we have 

a{2K + 2)a(l) = a{2K + l)a(2) + a{K + 2) + a{K + 1) 
a{2K + 2)-l = 1-1 + 1 + 1 
a{2K + 2) = 3 

Now we assume, as the inductive hypothesis, that a{2K + i) = —1 + 2i. We need to show 
that a{2K + {i + 1)) = -1 + 2(z + 1) = 2z + 1 for 2 < z < /sT. We know that 

a{2K +{i + l))a{i) = a{2K + i)a{i + 1) + a{K + {i + 1)) + a{K + i) 

Since 2 < i < K, we have 

3< i + 1 <K + 1 
K + 3< K + {i + l) <2K + 1 
K + 2< K + i <2K 

Therefore a{i + 1) = a{K + {i + 1)) = a{K + i) = l. 

a{2K + (i + 1)) ■ 1 = (-1 + 2i) ■ 1 + 1 + 1 
=^ a{2K + {i + I)) = 2i + l 
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So by induction, a{2K + i) = -1 + 2i for 2 < i < K + 1. 

Proofs of the other initial relations follow the same series of steps, and will therefore 
be left up to the reader to verify. □ 

Lemma 2. If the sequence {a{n)} is given by the quadratic recurrence then it also 
satisfies the linear recurrence given by 

a{n) = [2K^ + 8K + 4] {a{n - 2K) - a{n - AK)) + a{n - 6K) (5) 

for all n > 6K + 2, where the initial 6K + 1 values are taken to be the first 6K + 1 values 
of the quadratic recurrence 0). 

Proof. First note that by uniqueness, proving the converse (i.e.- that the sequence given 
by the linear recurrence ([5]) satisfies the quadratic recurrence (jll)) is equivalent to proving 
that the statement itself. Thus we prove that the linear recurrence given by ([5]) satisfies 
the quadratic recurrence given by (jlj) using strong induction. Define the sequence {a{n)} 
recursively for all n > 6K + 2 by ([5]) and let a{n) for 1 < n < 6K + 1 be defined by the 
initial conditions given in Lemma [H For simplicity in notation, let the term 2K'^ + 8K + A 
be called A{K). The proof will use methods analogous to those used by Hal Canary [2j 
to prove integrality of the Dana Scott recurrence ([1]). 
Let 

0(n) = a{n)a{n - {2K + 1)) - a{n - l)a{n - 2K) - a{n - K) - a{n - K - I) (6) 

We wish to prove by induction that 0(n) = for all n G Z+. Clearly for 1 < n < QK + 1, 
(t){n) = since the first 6K + 1 terms are defined to be the first terms given by (jlj). 

For the base case we must prove that (j){6K + 2) = 0. This is nothing but algebraic 
calculations easily verified by a computer algebra system such as Maple or Mathematica. 

0(6ir + 2) = a{6K + 2)a{4:K + 1) - a{6K + l)a{4:K + 2) - a{5K + 2) - a{5K + 1) 

where we can substitute the initial conditions for the linear recurrence for all but a{6K+2). 

a{AK + 1) = a(3Js: + {K + 1)) = 2K'^ + AK + I 

a{QK + 1) = a{5K + {K + 1)) = AK^ + 2AK^ + 40/^^ + IQK + 1 

a{AK + 2) = QK"^ + IQK + 9 

a{bK + 2) = AK'^ + 2AK'^ + A2K + 15 

a{bK + 1) = a{AK + {K + 1)) = AK^ + QK^ + lOiT + 1 

And a{QK + 2) must be given by the linear recurrence ([5]) 

a{QK + 2) = {2K'^ + 8K + A){a{AK + 2) - a{2K + 2)) + a(2) 
= {2K'^ + SK + 4)((6ir2 + IQK + 9) - 3) + 1 
= 12K^ + 80i^^ + 164^2 + 112K + 25 
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Now it is a matter of plugging this into your favorite computer algebra system and veri- 
fying that (p{6K + 2) = 0. 

Since the base case is verified we can proceed with the induction. We make the strong 
induction assumption that 0(m) = for m < n. We need to show that 0(n) = given the 
inductive hypothesis. Now, compute 0(n) by substituting for a(n), a{n — 1), a{n — K), 
and a{n — K — 1) from the definition of {a{n)}. 

a(n) = A{K)a{n - 2K) - A{K)a{n - 4K) + a{n - 6K) 
a{n - 1) = A{K)a{n -2K-1)- A{K)a{n - AK - I) + a(n - 6X - 1) 
a{n -K)^ A{K)a{n - 3X) - A{K)a{n - 5K) + a{n - 7K) 
a(n-K-l)^ A(K)a(n - 3K - 1) - A{K)a{n - 5K - 1) + a{n -IK -I) 

After substituting into expand and then simplify and we are left with 

4){n) = - A{K)a(n - 4:K)a{n -2K-l) + A{K)a{n - AK - l)a{n - 2K) 

- A{K)a{n - 3K) - A{K)a{n - 3K - 1) + a(n - 2K - l)a{n - 6K) 

- a(n - 2K)a{n - 6K - I) + A{K)a{n - 5K) + A{K)a{n - 5K - 1) 

- a{n - 7K) - a{n -7K -1) 

You can see that 

-A(K)a(n - AK)a(n - 2K - 1) + A{K)a{n - AK - l)a{n - 2K)+ 

- A(K)a(n - 3K) - A(K)a(n - 3X - 1) = -A(K)(f)(n - 2K) 

which equals by the induction hypothesis. Thus 

(t){n) = a{n -2K- l)a{n - 6K) - a{n - 2K)a{n - 6K - 1) + 

+ A{K)a{n - 5K) + A{K)a{n - 5K - 1) - a{n - 7K) + 
- a{n -7K-1) 

Substitute for a{n — 2K) and a{n — 2K — 1) from the definition of {a{n)}. 

a{n - 2K) = A{K)a{n - AK) - A{K)a{n - 6K) + a{n - 8K) 
a{n -2K-1) = A{K)a{n - AK - 1) - A{K)a{n - 6K - 1)+ 
+ a{n -8K-1) 

Simplify again to obtain 

0(n) = A(K)a(n - AK - l)a(n - 6K) - A(K)a(n - AK)a(n - 6K - 1)+ 

- A{K)a(n - 5K) - A{K)a{n - 5K - 1)+ 

+ a{n - 6K)a{n - 8K - 1) - a{n - 6K - l)a{n - 8K) + 

- a(n - 7K) - a{n -7K-1) 
= - A{K)(t){n - AK) + (t){n - Q>K) 
= 

Thus by induction 0(n) = for all n e Z+. □ 
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Theorem [T] now follows directly from the two lemmas. The first 6K + 1 terms of the 
sequence generated by the quadratic recurrence must be integers since they are given by 
polynomials with integer coefficients in the variables K and i. The rest of the terms in the 
sequence satisfy a linear recurrence with integer coefficients and integer initial conditions. 
Thus the terms must all be integers. 

4 Laurentness 

Now, instead of declaring that the initial conditions of a recurrence, r(n), be all 1, we will 
set r{i) to be the variable Xi for the first sufficiently many terms. If all terms in a sequence 
are Laurent polynomials in the initial variables {xi} then we say that the sequence (or 
the recurrence producing the sequence) has the Laurent property. In Section [H we gave 
examples of some sequences which were conjectured to be integers, and that have been 
shown to possess the Laurent property. Notice that a sequence having the Laurent prop- 
erty immediately proves integrality of that sequence. Since definition ([2]) is equivalent to 
defining the set of Laurent polynomials to be the subset of rational functions in which 
the denominator is a monomial, setting all variables equal to 1 in a Laurent polynomial 
clearly produces an integer. For this reason, and because the Laurent property has been 
helpful in producing a combinatorial proof of integrality in some cases, many people pre- 
fer proofs of Laurentness versus other proofs of integrality. In [3] , Fomin and Zelevinsky 
give easily verifiable sufficient conditions for a recurrence to possess the Laurent property. 
Using these conditions one can prove that the sequence given by (jl]) with initial conditions 
a{i) = Xi for 1 < i < 2K + 1 has the Laurent property. For the proof see pj. We do 
not include the proof here because the machinery used, namely Fomin and Zelevinsky's 
conditions via cluster algebras, are more advanced than our proof via linear recurrence. 
The technique of showing that a sequence generated by a quadratic recurrence also sat- 
isfies a linear recurrence, as in the proof from Section [3l has the potential to be applied 
in instances where the recurrence in question does not satisfy Fomin and Zelevinsky's 
conditions but still possesses the Laurent property, or even when the recurrence satisfies 
integrality without Laurentness. 

5 Generalization of this Family of Quadratic Recur- 
rences 

There is evidence to suggest that an even broader family of quadratic recurrence relations 
produces integer sequences. Let the sequence {b{n)}'^^^ be generated by the following 
recurrence. 

b{n)b{n — k) = h{n — i)b{n ~ k + i) + b{n — j) + b{n — k + j) (7) 

with the conditions that i<k — i<k, j<k— j<k, and b{l) = 1 for all / < k. With 
certain constraints on k,i, and j, the sequence is conjectured to be an integer 

sequence. 
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Conjecture 1. Consider the general form of the quadratic recurrence ^ with initial 
terms b{l) = 1 for < I < k. 

• In the case where k is even: 

— If i is odd, then j = | defines a recurrence that generates only integers. 

— If i is even, then j = f ; J = |; and j = ^ define recurrences that generate 
only integers. 

• In the case where k is odd: 

— If i is odd then j = ^ defines a recurrence that generates only integers. 

— If i is even then J = | defines a recurrence that generates only integers. 

Furthermore, all other values of j do not define a recurrence that gives integers exclusively. 

For example, the recurrence (jlj) is a special case of the generalization ([7]) where k is 
odd, i = 1, and j = in agreement with Conjecture [H 

6 Lifting the recurrences to 2-dimensional space 

As mentioned in Section |U these recurrences satisfy even stronger conditions than in- 
tegrality. However, the Laurent polynomials which are produced will have increasingly 
large coefficients. For the purposes of extracting combinatorics, it is preferable to modify 
these recurrences so that the sequences are Laurent polynomials with all coefficients of 
1 (which we will call faithful polynomials) so that each term counts some object (e.g. a 
perfect matching). 

Conjecture 2. For all conjectured families of recurrences in section 4 except for the case 
where k is even and j = |, the two dimensional recurrence 

Tin, k)T{n-k, k) = T{n-i, k + 2)T{n - k + i, k - 2) + 

+ T{n - J, A; + 1) + T{n -k + j, k - 1) 

with initial terms T{i,j) = Xij for all i < will generate faithful Laurent polynomials in 

Note that when k is even and j = |, we still obtain a Laurent polynomial with this 
method, but by combining linear terms one can see immediately that it is not faithful. 

Though it is not entirely clear what this lifting process will give us, it has been a 
means for finding combinatorics in similar sequences such as the Reid Barton and Dana 
Scott recurrences [1]. 
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7 Conclusion 



The integrality proof given here only sheds light on a small fraction of this new family 
of quadratic recurrences. It is our hope to eventually prove that all of the sequences 
generated by the quadratic recurrences in Conjecture [1] satisfy linear recurrences, and 
thus are integer sequences. The fact that the Laurent property is also satisfied in this 
special case (and possibly in the other cases of Conjectured]) may lead to the discovery of 
some family of combinatorial objects which are counted by this family of integer sequences. 

8 Acknowledgements 

We would like to thank James Propp for suggesting possible paths to take in this research 
and for running the Spatial Systems Laboratory at the University of Wisconsin-Madison 
where all of our research was conducted. Thanks also goes to the rest of the Spatial 
Systems Laboratory members for their helpful contributions along the way. We would 
also like to acknowledge and extend thanks to the NSF's Research Experiences for Un- 
dergraduates program and the NSA for funding our research. Finally we would like to 
thank Doron Zeilberger of Rutgers University for helping us make the final leap in the 
formulation and proof of Lemma [H 

References 

[1] Reid Barton. Personal communication, 2004. 
[2] Hal Canary. Personal communication, 2004. 

[3] Sergey Fomin and Andrew Zelevinsky. The Laurent phenomenon. Adv. in Appl. Math, 
28(2):119-144, 2002. 

[4] D. Gale. Mathematical entertainments: The strange and surprising saga of the Somos 
sequences. Mathematical Intelligencer, 13:40-42, 1991. 

[5] Emilie Hogan. An infinite family of quadratic recurrences satisfying the Laurent phe- 
nomenon. In preparation, 2008. 

[6] Michael Somos. Somos 6-sequence. http://grail.csuohio.edu/~somos/somos6.html. 

[7] Eric Weisstein. Laurent polynomial, http:/ /www. mathworld.com/LaurentPolynomial. html. 
From MathWorld-A Wolfram Web Resource. 

[8] Doron Zeilberger. Maple program "EmiliePaul(k,n)". 

http: / /www. math. rutgers.edu/~zeilberg/tokhniot/EMILIE. 



8 



